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Abstract. 

We present new mid-infrared spectra for a sample of 15 targets (1 FU Orionis object, 4 Herbig Ae stars, 5 T Tauri stars and 
O ' 5 Vega-type stars), obtained with the TIMMI2 camera at La Silla Observatory (ESO). Three targets are members of the yS Pic 

moving group (HD 155555, HD 181296 and HD 319139). PAH bands are observed towards the T Tauri star HD 34700 and 



i-i 



^ . the Herbig Ae star PDS 144 N. For HD 34700, the band profiles indicate processed PAHs. The spectrum of the Vega-type 

' object 7/ Corvi (HD 109085), for which a resolved disk at sub-mm wavelengths is known, is entirely stellar between 8-13 fim. 

Similarly, no indication for circumstellar matter at mid-infrared wavelengths is found towards the Vega-like stars HD 3003, 
HD 80951, HD 181296 and, surprisingly, the T Tauri system HD 155555. 
^ ' The silicate emission features of the remaining eight sources are modelled with a mixture of silicates of different grain sizes 

oo : and composition. Unprocessed dust dominates FU Ori, HD 143006 and CD-43 344. Large amorphous grains are the main 

' dust component around HD 190073, HD 319139, KK Oph and PDS 144 S. Both small grains and crystalline dust is found 

for the Vega-type HD 123356, with a dominance of small amorphous grains. We show that the infrared emission of the binary 
HD 123356 is dominated by its late-type secondary, but optical spectroscopy is still required to confirm the age of the system and 
' the spectral class of the companion. For most targets this is their first mid-infrared spectroscopic observation. We investigate 

, trends between stellar, disk and silicate properties and confirm correlations of previous studies. Several objects present an 

• exciting potential for follow-up high-resolution disk studies. 

^ , Key words. Circumstellar matter - Planetary systems: protoplanetary disks - Stars: pre-main sequence - Infrared: stars 

d 1. Introduction tures of 27 T Tauri stars (TTS) and found that the overall degree 

of crystallinity in T Tauri disks correlates neither with stellar 
Up to now a full understanding of dust evolution in pre-main age nor luminosity. Through resolved interferometric observa- 
sequence (PMS) disks and the survival of debris disks is lack- tions, van Boekel et al. ((2004J showed that a higher amount of 
ing. Both disk types need to be distinguished in terms of their crystalline dust is seen in close proximity to the star than in the 
dust evolution. For an in-depth discussion of gas and dust evo- outer disk regions. In an anal ysis of 40 TTS and 7 HAeBe stars, 
lution in circumstellar disks, we refer to the review by Meyer Kessler-Silacci et al. (120061) found a correlation with spectral 
et al. (I2007I I. type for the ratio strength-to-shape in the N-band silicate fea- 

ture, but no relationship with stellar age or disk evolutionary 
Primordial PMS disks largely reflect the grain size distribution ^^^^^^ Therefore, the observed grain sizes in PMS disks must 

of the interstellar medium (ISM), while more evolved PMS i j j »u * . i-i i 

^ ' also depend on other factors luce, e.g., close companions, or 

disks experience grain growth and crystallisation. However, t^^^buience and gas content in the disk, while the exact depen- 

there is no global correlation of dust properties with age, not ^^^^.^^ ^^^^^^^ ^.^^^^ unknown. Kessler-Silacci et al. (l2007l 

evenwhen considering individual spectral types. Schegerer et demonstrated that the disk region which is probed by the 10 /zm 

al. re-analysed th e mid-infrared silicate emission fea- ^-^-^^^^ ^^^^^i^j^ ^^^^^^^^^ ^^1^^^^ ^^^U^ luminosity, what 

Send offprint requests to: oschuetz@eso.org '^o^ld explain the correlation of grain size with spectral type. 

* Based on observations made at the European Southern Obser- 
vatory, obtained under program IDs 076.C-0634(A) and 077.C- Debris disks - also referred to as Vega-type objects - are of- 
0054(A), during TIMMI2 technical time as well as collected from the ten characterised by the lack of gas, so that the dust dynam- 
ESO/ST-ECF Science Archive Facility ics is no longer influenced by that of the gas. The presence 



Table 1. Stellar parameters and known fluxes of our target sample. 



Object 


Class 


Spectral Type 


T 1 , H 

A adopted 


Hop" P"^ 1 „f J 

V^'-'fe 5 ^adopted 


V 


Fn 

rjvi 


Ape 
rMvrl 


Ref. 


d 


Ref. 


A,, 
rtnas?! 


Ref 


Sep. 

L J 


Ref. 


HD 3003 


Vega 


AOV 


11000 


4.5 


5.07 


0.48 


50 


(1) 


47 


(9) 






0.1 (?) 


(26) 


HD 34700 A 


TTS 


GOIVe 


6000 


4.5 


9.15 


0.60 








(10) 


0.68 


(19) 


SB/5.2/9.2 


(27) 


HD 80951 A 


Vega 


AlV 


9500 


4.5 


6.00 


0.44 






220 


(9) 






0.3/7.1/48.2 


(26) 


HD 109085 


Vega 


F2V 


7000 


4.5 


4.31 


2.18 


1000 


(2) 


18 


(11) 


0.12 


(20) 






HD 123356 AB 


Vega 


GIV/M? 


6000/4000 


4.5 


10.0/12.2 


1.36 








(12) 


0.16 


(20) 


2.5 


(28) 


HD 143006 


TTS 


G5Ve 


5750 


4.5 


10.21 


0.86 


~1 


(3) 


82 


(13) 


0.62 


(21) 






HD 155555 ABC 


TTS 


G5IV/K0IV-V/M4.5 


5750/... 


4.5 


6.88/13.0 


0.69 


12+* 


(4) 


32 


(9) 






SB/ 33.0 


(26) 


HD 181296 A 


Vega 


AOVn 


9500 


4.5 


5.03 


0.54 


-4 


(4) 


48 


(9) 


0.12 


(20) 


4.1 


(29) 


HD 190073 


HAeBe 


A2IIIpe 


9500 


4.5 


7.82 


7.16 


1 9 +0.6 
^•^ -0.6 


(5) 


>290 


(14) 


0.19 


(19) 






HD 319139 


TTS 


K5Ve / K7Ve 


4250 


4.5 


10.44 


0.45 


5 5 +'■' 

J.J J 


(6) 


83 


(15) 


< 0.12 


(22) 


SB 


(30) 


CD-43 344 


TTS (?) 


M2 


4500 


4.5 


9.42 


3.73 








(10) 


0.0 


(20) 






FU Ori A 


FUOR 


Gon 




4.5 


8.94 


5.95 


<0.3 


(7) 


450 


(16) 


2.2 


(23) 


0.5 


(31) 


KK Oph A 


HAeBe 


A6Ve 


8750 


4.5 


9.4-12.9 


9.87 


7 +1.0 
' -0.5 


(8) 


160 


(17) 


1.6 


(24) 


1.6 


(32) 


PDS 144 N 


HAeBe 


A2IV 


8750 


4.5 


14.2 








-1000 


(18) 


3.0 


(25) 


-5.0 


(33) 


PDS 144 S 


HAeBe 


A5V 


8750 


4.5 


13.1 








-1000 


(18) 


1.1 


(25) 


-5.0 


(33) 



References: Spectral types are obtained either from SIMBAD or one of the cited papers, if a more accurate classification is available in the literature. Tadoptcd and (log g)adoptcd list the temperature 
and surface gravity which we used for the Kurucz atmosphere models. V-band magnitudes and IRAS 12 fim fluxes are taken from the SIMBAD database. KK Oph shows variability in the 
visual. PDS 144 is not resolved by IRAS; its V-band fluxes were given by Perrin et al. ( I2006t . Stellar ages: (1) Song et al. (2001^), (2) Wyatt et al. (2 005, 1 and references therein, (3) Dunkin et 
al. l |1997t and references therein, (4) mean age of the yS Pic moving group as given in Zuckerman & Song (2001), (5) Catala et al. ( 12007 1 . (6) Quast et al. ( I2000I I, (7) e.g. D'Angelo et al. ( 120001 , 
(8) Carmona et al. ( 120071 . Errors in the age determinations are quite heterogeneous, depending on the cited work. Distances are taken from (9) Hipparcos parallax, (10) uncertain Hipparcos 
parallax, (11) Wyatt et al. ( 120051 , (12) The Hipparcos parallax was recently removed from SIMBAD and is likely not firm, (13) Sylvester et al. ( 119961 , (14) van den Ancker et al. ( 119981 , 
(15) Quast et al. (2000), (16) e.g. Wang et al. (2004), (17) Hillenbrand et al. ( 11992 1 (18) Perrin et al. (2006 1 and references therein. Extinctions: (19) van den Ancker et al. (1998), (20) Sylvester 
& Mannings ( 2000), (21) Malfait et al. {199S), (22) Hutchinson et al. ( fT990t , (23) Kenyon et al. (JJ88 I, (24) Carmona et al. (2007), (25) Perrin et al. ( |2006l l. Stellar separations in multiple 
systems: (26) Dommanget & Nys i2002t . (27) Sterzik et al. i2005] |, (28) Worley & Douglass ([1997}, (29) Lowrance et al. ( 2000), (30) Quast et al. ( [20001 , (31) Wang et al. ( [20041 . (32) Leinert 
et al. ( 119971 , (33) Perrin et al. (1 20061 . 'SB' denotes a spectroscopic binary. Note on the multiplicity of HD 155555: Although the SB components are not spatially resolved, we label them as A 
and B, since this naming convention existed already in the literature. 
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Object 


Origin 


Airmass (N-spec) 


tin, (N-spec) 


Fni 


Fn10.4 


Fni 1.9 


Ref. 


N-spec features 








[min] 


[Jy] 


[Jy] 


[Jy] 






HD 3003 


technical 


1.2 


24 


0.58 ± 0.03 


0.32 ± 0.03 


0.20 + 0.10 




stellar 


HD 34700 


technical 


1.2-1.5 


60 


... 


... 


0.79 ± 0.29 




PAH 


HD 80951 


077.C-0054 


1.7 


53 


... 




0.44 ± 0.03 


(1) 


stellar 




cllt-Ill VC 


1.1 


12 






1.51 + 0.24 




t;tf*l 1 Qi" 
SLCilcli 


HD 123356 


archive 


1.0 


18 




1.51 ±0.05 


1.68 ±0.06 




silicate emission 


HD 143006 


archive 


1.1 


24 






0.86 ± 0.06 


(1) 


silicate emission 


HD 155555 


077.C-0054 


1.3 


58 






0.46 ± 0.03 




stellar 


HD 181296 


077.C-0054 


1.2 


64 


0.59 ± 0.05 








stellar 


HD 190073 


technical 


1.5 


12 






7.16 ±0.43 


(1) 


silicate emission 


HD 319139 


077.C-0054 


1.0 


160 






0.36 ± 0.03 




silicate emission 


CD-43 344 


technical 


1.1 


24 


3.78 ± 0.05 


3.82 ± 0.05 


3.04 ± 0.05 




silicate emission 


FUOri 


technical 


1.4 


20 




6.00 ± 0.03 


5.61 ± 0.03 




silicate emission 


KKOph 


archive 


1.2 


12 


10.09 + 0.05 




10.71 ± 0.08 




silicate emission 


PDS 144 N 


077.C-0054 


1.0 


32 






1.74 ± 0.09 


(3) 


PAH 


PDS 144 S 


077.C-0054 


1.0 


13 






2.49 ±0.13 


(3) 


silicate emission 



Results from TIMMI2 photometry and spectroscopy are merged in this table. The airmass and TIMMI2 integration times refer to N-band 
spectroscopy. Mid-IR photometry was obtained in the specified TIMMI2 filters with central wavelengths 8.6 /jm (Nl), 10.3 /jm (N10.4) and 
1 1.6 yum (Nl 1.9). Errors represent the accuracy of the aperture photometry. Where no own measurements were available, e.g, for some archive 
data, we adapted data from (1) IRAS, (2) Wyatt et al. ( 12005b or (3) Perrin et al. ( 12006b . 



of small dust particles in such disks points towards a contin- 
uous replenishment of the dust grains, e.g. by collisions be- 
tween planetesimals, asteroids or comets. Indeed, through col- 
lisional cascades significant amounts of small grains can be 
produced, down to grain sizes which are easily removed by 
radiation pressure on short timescales, and therefore, without 
replenishment would not be expected to be present in the disk. 
Song et al. (2005) found a surprisingly strong silicate feature 
for the -300 Myr old star BD-i-20 307, which they explain by 
the release of secondary dust after the destruction of asteroids. 
The emission peak shortward of 10 jum in its spectrum indi- 
cates rather small grains. In Schiitz et al. (|2005b hereafter 
Paper II), we discussed a similar - but much younger - main- 
sequence object, HD 113766, which shows large amounts of 
crystalline secondary dust, possibly released by collisions in a 
planetesimal belt. By applying spatially resolved spectroscopy 
to jS Pic, Okamoto et al. j2004l) were able to trace the origin 
of secondary dust towards the location of putative planetesi- 
mal belts. HD 69830 is another system with secondary debris 
dust in form of warm, small, crystalline silicates (Beichman et 
al. 12005b . possibly resulting from the disruption of an aster- 
oid within 1 AU from the star (Lisse et al. 2008). Interestingly, 
this system is also known to host three Neptune-sized planets 
(Lovis et al. l2006] l. 



The present paper continues our previous mid-infrared (MIR) 
spectroscopic studies of circumstellar disks. In Schiitz et al. 
(I2005al hereafter Paper I) we described our strategy and dust 
emission models, and analysed a sample of pre-main sequence 
stars. Vega-type objects with debris disks were studied with the 
same methods in Schiitz et al. (I2005bl Paper II). 



This contribution describes new mid-infrared spectra of cir- 
cumstellar disks. We selected targets with infrared excess that 
have not yet been discussed in depth in the literature, i.e. ob- 
jects which were suspected to be disk candidates, but for which 
no MIR spectroscopy had been obtained before. Several iso- 
lated targets were unobserved by Spitzer, as they were not in- 
cluded in the surveys of open clusters and moving groups. We 
re-observed targets where newer or higher SNR data were re- 
quired, as for HD 34700. A MIR spectrum of KK Oph had al- 
ready been published, however, it was not analysed with the 
methods we present here. We included the well studied object 
FU Ori to enable a comparison with other FU Orionis targets 
from Paper I. Table [1] summarises our targets and their stellar 
parameters. In Sect.|2]we describe the observations and data 
reduction. Sect. [3] discusses the spectral energy distributions 
(SEDs) and methods for analysing the emission feature in the 
8-13 jum spectra. In Sect.|4]we describe observed trends in the 
dust properties. A summary of our results is given in Sect.|5] 
The known properties of each target are described in detail in 
the appendix. It is our goal to contribute a piece to the puzzle 
of disk and dust evolution around young stars. 



2. Observations and data reduction 

TIMMI2: The principal data described in this work were ob- 
tained with the ESO TIMMI2 camera (Kaufl et al. 12003]) at 
La Silla Observatory. We performed most of the observations 
during run 077.C-0054(A) in June 2006, supported by fur- 
ther data retrieved from the ESO/ST-ECF Science Archival 

' http;//archive.eso.org 
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0.1 1.0 10.0 100.0 1000.0 



X [|im] 

Fig. 1. Spectral energy distributions for FU Ori and the 
Herbig Ae stars. A Kurucz model is overplotted for all tar- 
gets to represent the stellar contribution to the SED, apart from 
FU Ori as justified in the text. The profile of the N-band spec- 
trum is indicated in the SEDs (at the position of the arrow) to 
show their agreement with the photometry. For PDS 144 N and 
S all IR fluxes are taken from Perrin et al. (I2006I I. 



which originally had been observed in March 2003 (tar- 
gets HD 109085, HD 123356 and HD 143006; ESO program 
070.C-0743(A); PI J.-C. Augereau) and July 2003 (target 
KKOph; ESO program 071.C-0438(A); PI H. Linz). Further 
data were collected during TIMMI2 technical nights in October 
2005 and January 2006. Table [U summarises the origin of all 
data. 

Imaging data were obtained with a 0.20"/pixel scale, corre- 
sponding to a 64" x48" field of view. For most targets the nod- 
ding direction was perpendicular to the chop throw, resulting 
in four images of the source in the co-added frame, each carry- 
ing 25% of the flux. Some archive data, however, had the nod 
parallel to the chop, giving three images of the source, with the 
central one carrying 50% of the flux. We applied aperture pho- 
tometry with radii rap between 2.6 pixel and 4.7 pixel, depend- 
ing on the target. The optimal radii were derived from curves 
of growth and averaged over the radii obtained for each of the 
beams, to account for detector inhomogeneities. As error in the 
photometry we give the mean of the measurement errors for 
each beam. The background level was determined by measur- 




10'- /CD -43 344 

10-'° U . . . . ^ 

0.1 1.0 10.0 100.0 1000.0 10000.0 

X [(im] 



Fig. 2. The spectral energy distributions for the T Tauri stars. 
Remarkably, HD 155555 shows no excess emission in the in- 
frared. 



ing the mean sky value per pixel around each beam in an an- 
nulus centered on the star with inner radius rap + 10 pixel and 
outer radius rap + 20 pixel (resulting in a maximum outer radius 
of 5", while the chopping amplitude is 10"). Photometric stan- 
dard stars were selected from a list of MIR standards by Cohen 
(119981) and observed approximately each 2 hours. 

N-band spectra were obtained between 8-13 yum applying a 
standard chopping and nodding technique along the 1.2" slit 
(for archive data and technical time) or the 3" slit (run 077. C- 
0054(A)) with a throw of 10". The spectral resolution of the 
used 10 jum low-resolution grism is 0.02 yum/ pixel {A ISA ~ 
160). Airmass and on-source integration times of the observa- 
tions are listed in Table[T] Standard stars for telluric correction 
and flux calibration were selected from the list of Cohen (1 19981 1 
and observed close in time and airmass. 

The spectroscopic frames were de-biased and checked for a 
correct slit alignment along the detector, i.e. whether the spec- 
trum is aligned with the x-axis. After co-adding the positive and 
negative nodding positions, the extraction profile is determined 
from the cross dispersion profile: the target spectrum is ex- 
tracted along a profile width where the SNR exceeds 2.5 times 
the noise of the background. We furthermore applied an air- 
mass correction described in Schiitz & Sterzik (I2005cl) and cor- 
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0.1 1.0 10.0 100.0 1000.0 

X [|im] 

Fig. 3. Spectral energy distributions of the Vega-type stars. For 
the binary HD 123356 the SED of the primary {black curve) is 
entirely stellar, while the excess emission is dominated by the 
secondary {grey curve and stars). 



relate the spectrophotometric calibration with the obtained N- 
band photometry. The wavelength calibration is optimised us- 
ing the atmospheric absorption features of ozone (9.58 //m) and 
CO2 (11.73 //m, 12.55 //m) in addition to the given TIMMI2 
wavelength calibration table, i.e. shifts were applied to min- 
imise residual division between the atmospheric features in 
standard stars and target spectra. When atmospheric conditions 
were unstable between observations of target and standard star, 
leaving uncorrectable ozone features around 9.5 /im, we cut 
these remnant features from the spectra. By comparing data 
from different nights and applying several standard stars for 
comparison, the exact position of the ozone remnant feature 
can be identified. Depending on the atmospheric quality, this 
leaves gaps of different sizes up to a maximum range between 
9.4 and 9.9 /vm. The SNR in our target spectra, averaged be- 
tween 8-12 jum, ranges between 10 and about 70 for all ob- 
jects, with a minimum SNR in the spectra of HD 3003 and 
HD 319139, while the highest SNR is obtained for FU Ori and 
KK Oph. 

Spitzer: HD 34700 is also observed with the IRS on board 
the Spitzer Space Telescope (Werner et al. |2004| l (AOR key = 
3584768, GTO program of J. Houck). The spectra were taken 
with the Short-Low (SL, 5-15 /im) and Short-High (SH, 10- 
20 jum) modules (Houck et al. I2004l l with a resolving power 



of respectively R = 60-128 and 600. The basic data reduc- 
tion was performed by the standard IRS pipeline S13 at the 
Spitzer Science Center (SSC). Subsequently, bad pixels were 
removed with IRSCLEAN (available from the SSC website|l) 
and a standard point-source extraction has been performed us- 
ing SMART (Higdon et al. |2004| l: tapered column extraction 
for SL data (i.e. the width of the extraction scales with wave- 
length and equals 7.2" at 6 //m); full aperture extraction for SH 
data (4.7" X 11.3"). 

Further auxiliary data: During the discussion of two targets 
- HD 109085 and HD 123356 - we will refer to individual 
data obtained either with VISIR (Lagage et al. l2004l l or SOFI 
(Moorwood et al. |1998l l. 



3. Analysis 

3.1. Spectral Energy Distributions 

Together with fluxes from the literature in the passbands 
UBVRI, JHK (2MASS), 12, 25, 60 and 100 fim (IRAS) we 
used the TIMMI2 photometry in Table [T] to construct a spec- 
tral energy distribution (SED) for our targets. In Figs. [TJ{3] the 
resulting optical to mm spectral energy distributions are dis- 
played. To emphasise the non-stellar contribution to the SED, 
we plotted a Kurucz atmosphere model (Kurucz 119941 ) with 
the stellar parameters from Table [T] normalised to the V-band 
flux. For FU Ori we did not plot a Kurucz model, as for this 
kind of object the V-band flux most likely arises from the star 
plus the hot inner part of the outbursting disk (cf. Sect. |A. It . 
Where required, we reddened the Kurucz model with the ex- 
tinction Av given in Table [T] The photometry to construct the 
SEDs was taken from Hutchinson et al. ( 1 19901 1, Henning et al. 
( 1994), Sylvester et al. (fT996b . Coulson et al. (fT998b . Malfait et 
al. (1998 ), Mannings & Barlow (119981) . Sylvester & Mannings 
(2000), Sylvester et al. (|200TT) . Jayawardhana et al. (|200TT ). 
Mamajek et al. (200?), Natta et al. ("2004), Sheret et al. (12004] ). 
Wyatt et al. (2005), Perrin et al. (i2006]) . the catalogs JPll, 
2MASS, IRAS, and SIMBAD. 

The SED of HD 155555 can entirely be explained by stellar 
emission, while for HD 3003, HD 80951, HD 109085 and 
HD 181296 the IR excess begins longward of 25 yum. We, 
therefore, do not consider those objects in our dust analysis. For 
the other targets a near-infrared (NIR) excess is seen, starting 
longward of the L-band for HD 34700 as well as HD 319139, 
and around 1 jim for the remaining eight targets. 

For HD 109085 (alias rj Crv) we obtained Ql band (17.7 jum) 
imaging with VISIR during run 076.C-0634(A) in March 2006. 
77 Crv was observed using a standard chopping technique with 
a throw of 10" at an airmass of 1.03, while the standard star 
(HD 99167) was obtained at airmass 1.35. Standard data re- 
duction and photometric extraction was applied, including a 
correction for the atmospheric extinction. The derived photom- 
etry is 1.59 Jy (+ 0.10 Jy). Subtraction of the PSF showed that 

- http;//ssc. spitzer.caltech.edu 
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Table 3. Relative mass fractions in percent derived from the model fits to the N-band spectra. "Amorphous silicates" indicate the 
sum of amorphous olivines and pyroxenes. Targets which do not show a silicate feature are not listed. 



Object 


Amorphous silicate 


Forsterite 




Silica 


Enstatite 




0.1 yum 


2.0 yum 


0.1 yum 


2.0 yum 


0.1 yum 


2.0 yum 


0.1 yum 


2.0 yum 


HD 123356 


76.0!|-3 


0+" - 

"•"-0.0 


' •"-o.s 




1 1 4^"^ 

11.1- Q 3 




^•"-0.6 


1 o+'-° 

^•"-1.0 


HD 143006 


76 4+^ '' 




3 0+"' 

-'■"-0.2 


1 9+"-3 




'■-'-0.7 


1 0+°^ 


8.2+0;^ 


HD 190073 


-'•'-0.6 


85.3^Q 2 








7 9+0.1 
'■"^-0.1 


T Q+0.1 

-■•'-0.1 




HD 319139 


9c 7+18.4 
^-'•'-17.0 


63.6^^ 


^^■^-0.4 


A 7+1-6 
^■^-1.7 


'■-'-i.i 


3.8!|:i 






CD-43 344 


74.3!»:J 






"■^-0.3 










FUOri 


87 9+^-^ 




"■-'-0.0 


'■"-0.1 


1 4+0.2 
^•^-0.1 




1 4+0.2 
^•^-0.1 




KKOph 


q 7+0.8 
^•'-0.8 




^■^-0.1 


"■^-0.4 










PDS 144 S 


30 7+'-'' 




"■-'-0.1 






14.2!«j 







the source is not resolved at 17.7 /im: the FWHM is 6.7 pixel 
(0.50") for the standard star while only 6.8 pixel (0.51") for 
rj Crv. 



We determined from which component of HD 123356 the ex- 
cess emission arises. The 2MASS photometry gives the com- 
bined fluxes of both sources but interestingly, the correspond- 
ing 2MASS coordinates seem to be centered towards the sec- 
ondary. With SOFI images obtained during thick clouds in 
May 2006 we could separate the emission and find that the 
secondary dominates the IR. The relative NIR brightness dif- 
ference between secondary and primary amounts to 3.3 mag 
in filter NB1.19, 4.0mag in NB1.71 and Fell-H, and 4.2 mag 
in NB2.09. With this information it was possible to split the 
2MASS fluxes between both components and to plot a resolved 
SED (see Fig.O. From the VJHK colors of both stars, the pri- 
mary appears as a F7-G1, in agreement with the SIMBAD clas- 
sification, while the secondary seems to be either a very late 
M-type or it sufl'ers extinction. In TIMMI2 MIR data we only 
see one isolated source. Although astrometry is not possible on 
TIMMI2, with the observations discussed above, we can safely 
assume that we only see the secondary in the MIR. In the SED 
we plot Kurucz models for both components. 



3.2. Dust properties 



In Figs.|4}{6] and [8] we show the N-band spectra for our tar- 
gets. To quantify the differences between the sources and to 
determine the composition of the circumstellar dust, we adopt 
the same procedure as described in Sect. 4 of Paper I, which 
is summarised below in Sect. |3.2.T] For further details we re- 
fer to Paper I. Since two objects show strong PAH emission, a 
description of PAH classes is given in Sect. l3.2.2l 



3.2.1. Silicates 

When stars form out of their parental cloud, the material in 
their disks is assumed to have a similar initial composition as 
in the interstellar medium, in which amorphous silicates are the 
main component observed at 10 yum (e.g. Kemper et al. l2004l l. 
Changes in the composition, structure and size are expected to 
occur during the subsequent evolution of the star-i-disk system, 
possibly leading to a planetary system. Laboratory experiments 
have suggested that, due to thermal annealing, amorphous sil- 
icates gradually turn into crystalline forsterite and silica (e.g. 
RietmeijerT989^; Hallenbeck & NuthTW). Bouwman et al. 
( 200 1 ) found a correlation between the amount of forsterite and 
silica in the disks of Herbig Ae/Be (HAeBe) stars, indicating 
that thermal annealing indeed takes place in these objects. A 
similar correlation was also found for the lower-mass T Tauri 
stars (e.g. Meeus et al. l2003l) . These authors also concluded that 
the dust around HAeBe stars and T Tauri stars has very similar 
characteristics. In Paper II, we successfully modelled a sample 
of Vega-type stars with the same dust species as those found in 
young stellar objects. The above-mentioned dust species emit 
in the N band, making this an excellent window to study dust 
evolution in the inner parts of the circumstellar disk. 

Our analysis is based on a method that is commonly used when 
analysing the dust features around 10 fim: fitting a continuum 
to the spectral range outside the feature, and then decomposing 
the continuum-subtracted emission feature by fitting a linear 
combination of absorption coefficients of different dust species, 
for a few selected sizes. For a deeper discussion on this method, 
we refer to Bouwman et al. (1200 11 1 and van Boekel et al. (I2005I I. 
To derive the composition of the circumstellar dust, we model 
the spectra in Figs. |4}{5] with a linear combination of emission 
features from the following dust species, which have been iden- 
tified in the disks of pre-main sequence stars: 

- Amorphous olivine ([Mg,Fe]2Si04) with grain sizes of 0.1 
and 2.0 fj.m, to which we will refer to as "small" and 
"large" grains. We used absorption coefficients of spheri- 
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Fig. 4. Dust decomposition for FU Ori and the Herbig Ae stars. 
The different linestyles - or colours in the online version - 
represent small amorphous silicates (dash-triple-dots, orange), 
large amorphous silicates (dashed, red), both small and large 
grains of crystalline forsterite (dash-dots, blue), silica (SiOa, 
dotted, purple) and enstatite (long curved dashes, green). The 
very thick black curve corresponds to the observed spectrum. 
The sum of the different model components is given by the 
overplotted grey curve. 



cal grains based on laboratory experiments by Dorschner et 
al. (fT995T l. 

Amorphous pyroxene ([Mg,Fe]Si03), as well with grain 
sizes of 0.1 and 2.0 fim. These absorption coefficients are 
also based on Dorschner et al. (|19951 l. 



Fig. 5. Dust decomposition for the T Tauri stars (first three pan- 
els). See Fig.|4]for an explanation of the linestyles. The bottom 
panel shows HD 123356, the only Vega-type object in our sam- 
ple that shows a silicate feature. 



Crystalline silicates: magnesium forsterite (Mg2Si04) and 
enstatite (MgSiOs). The latter one, however, is a rare 
find even in young sources (cf. Paper I). Coefficients for 
forsterite are from Servoin & Piriou (I1973I I. for enstatite 
from Jaeger et al. ( |19981 l. Each of these silicates was mod- 
elled with "small" and "large" grains. 

Silica (SiOa) coefficients were taken from Spitzer & 
Kleinman ( 11961) . 
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Shape and size of the dust grams: Astronomical dust grains 
likely have various shapes and structures (compact/fractal), 
which are reflected in their different spectral signatures, com- 
plicating the spectral analysis. Min et al. (2003, 2005) showed 
that mass absorption coefficients, calculated with a distribution 
of hollow spheres (DHS), give a reliable identification of solid 
state components. In this work we use mass absorption coeffi- 
cients Ki from van Boekel et al. ( 120051) . More specifically, for 
the crystalline grains and silica, a-, was calculated with DHS, 
while for the amorphous olivine and pyroxene grains Mie the- 
ory was applied, as here shape effects are less significant. For 
a more extensive discussion of shape and structure effects on 
spectral signatures, we refer to Min et al. (2003, 2005). 

In circumstellar disks, dust grains grow and get destroyed by 
various processes. As a result, the dust grains do not have a 
single size, but a wide range in sizes. Since the grains are not 
always compact, but can be fluffy as well, the term 'volume 
equivalent radius', ry is introduced: it is the radius of a sphere 
with the same material volume as the particle. This means that 
fluffy particles have a much smaller ry than their real size. Van 
Boekel et al. (12005b showed that the size distribution of grains 
radiating in the 10 fim region can be represented by two grain 
sizes: 'small' (ry = 0.1 jum) and 'large' {ry - 1.5 yum). It is this 
approach that we follow when modelling our spectra. 

Fitting procedure: We used absorption coefficients from van 
Boekel et al. ( I2005I I and Bouwman ( 120081 1. which are based 
on the laboratory data cited above. For the compositional fit 
we perform an error analysis as described in Sect. 5.1.3 of van 
Boekel et al. (|2005] |. by minimising the expression for the re- 
duced 

1 ^-1 (Trmodeli' i \ rTrobserved/ i \ ^ 
^.ed - h 

As contains 12 free parameters, i.e. the mass fractions 

of the five silicate emissivities for two grain sizes each and the 
continuum and dust temperature, M equals 12. The continuum 
temperature was determined by fitting a blackbody to the con- 
tinuum. A^^ is the number of wavelength points /I,, while cr, 
describes the absolute error of the observed flux at wavelength 
A, and corresponds to the background noise derived from the 
data. The eiTors of the silicate mass fractions are derived from 
a Monte-Carlo simulation. For each spectrum we create 100 
synthetic spectra by randomly adding Gaussian noise with a 
width up to the level of the background fluctuations at each 
wavelength point. The same compositional fitting procedure is 
applied to each of these spectra, resulting in slightly different 
mass fractions, from which the mean values and errors are ob- 
tained. 

We caution against overinterpreting our fitting results. Indeed, 
significant systematic errors could be introduced by either an 
imperfect calibration or by the limited wavelength range and 
SNR, both inherent to ground-based observations (see Juhasz 
et al. 120091 for a thorough discussion on this topic). Since we 
cannot quantify these systematic errors, we could also not take 
them into account in the fitting procedure. Nevertheless, the re- 
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Fig. 6. The TIMMI2 spectrum of PDS 144 N, showing the 8.6, 
1 1 .0 and 1 1 .2 nm PAH bands. 



suits are robust in terms of distinguishing between small and 
large grains, but probably only indicative in terms of the crys- 
talline species. 

The modelling results for the emission spectra, together 
with the contribution of each dust component are shown in 
Figs-EHl] Note that, while we model the spectra with the 10 
(grain and size) dust types and PAH as explained above, we 
only plot 6 components for a better visibility. In particular, we 
added amorphous olivine and pyroxene of 0. 1 ^va and plot them 
as "small amorphous silicates". Similarly, "large amorphous 
silicates" contain the 2.0 yum pyroxene and 2.0 yum olivine 
grains. For silica, forsterite and enstatite, we plot the sum of 
the 0. 1 and 2.0 /zm grains. 

When deriving in what amount the different dust species are 
present, the linear coefficients of the fit - which are propor- 
tional to the radiating surface of the grains, assuming optically 
thin emission - need to be converted to mass. However, it is 
not possible to determine the absolute amount of mass for each 
species present, as we have no spatially resolved data to de- 
rive the particles' size, density or temperature distribution as a 
function of radius within the disk. Therefore, we derived mass 
fractions of the different species under the assumption that (1) 
the particles are spherical and (2) the particles have the same 
density. The mass fractions are meaningful to compare the ob- 
jects in our sample and to establish the relative amount of pro- 
cessed dust we observe at 10 im\. In Table [3] we list the de- 
rived mass fractions, which is the ratio of the mass of a partic- 
ular dust species and the total dust mass. Please note that for a 
given mass, small particles have a larger total emitting surface 
than large ones. Inversely, a similar amount of observed radia- 
tion will result in a much smaller mass when caused by small 
grains than if it was caused by larger grains. It is also impor- 
tant to point out that our results are only valid for the warm 
dust, which is located in the inner part of the disk atmosphere 
(r < 10 AU) and radiates in the 10 yum region. If the disk atmo- 
sphere and midplane are well-mixed, the relative abundances 
would reveal the bulk composition of the inner disk, but if there 
is vertical settling the interpretation is more complex. 
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Fig. 7. A comparison of the 8.6 and 1 1 .2 //m PAH band profiles 
of PDS 144 N (black) with the PAH classes (grey shades) as 
determined by Peelers et al. (2002) and van Diedenhoven et al. 
(I2004i) . Clearly, PDS 144 N has a "class A" 8.6 ^im PAH band, 
while the 11.2 fim PAH profile is different from both class A 
and B. The wavelength error is 6 A/ A ~ 10-^ for the science 
spectra, and 6 A/ A ~ 10 for the PAH class model profiles. 



3.2.2. PAH 
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Fig. 8. TIMMI2 (grey) and Spitzer-IRS (thin black and thick 
black) observations of HD 34700 are shown. PAH emis- 
sion bands are clearly present on top of the dust continuum. 
Atmospheric variations may have influenced the spectrophoto- 
metric accuracy of the TIMMI2 data. 



HAeBe stars and Post-AGB stars reveals a correlation of the 
position of the 7.7 and 11.2 /urn PAH bands with the effective 
temperature of the central star (Sloan et al. 120051 120071) . 



Polycyclic Aromatic Hydrocarbons (PAHs) show N-band 
emission features at 7.7, 8.6, 11.0, 11.2 and 12.7 fim. A key 
result in observational studies of the PAH bands is that their 
profiles show clear variations in terms of peak positions, pro- 
file shape and relative intensities (Peeters et al. |2004al and ref- 
erences therein). HII regions, reflection nebulae and the ISM 
all show the same PAH band profiles (called "class A"), while 
planetary nebulae and post-AGB stars have clearly different, 
red-shifted profiles, most noticeable for the 6.2 and 7.7 yum 
emission bands (class B/C). Figs.|2]and|9]show these PAH pro- 
files for two of our targets. Note that in contrast to class A, 
in which the PAH profiles in all objects are almost identical, 
class B comprises a large variety of profiles for the 6.2 and 
7.7 fim PAH bands (Peeters et al. 120021 ). The latter class is 
based on the spectra of evolved stars, where PAH formation 
and processing is actively taking place. Also HAeBe stars show 
large variations in their PAH profiles (Peeters et al. 2002; van 
Diedenhoven et al. 2004; Sloan et al. ,2005) . HAeBes with pro- 
files similar to those found in the ISM (class A) turn out to 
be embedded in their parental cloud and PAHs from this cloud 
are dominating the emission. Thus, the initial characteristics of 
the PAH family entering protoplanetary disks of HAeBe stars 
are the same as in the ISM. On the other hand, those HAeBes 
exhibiting significantly different PAH profiles (class B) are iso- 
lated HAeBe stars and spatial studies show that the PAH emis- 
sion originates in the disk (e.g. Habart et al. 120041 120061 Geers 
et al. I2007bl Lagage et al. i2006t Doucet et al. |2007] |. imply- 
ing that PAH processing is actively occurring in these disks. 
In addition, the presence of PAH emission is related with the 
disk geometry (Meeus et al. 1200 U Acke & van den Ancker 
|2004| l: PAH emission is observed in flaring disks where the 
PAH molecules can be excited by UV photons from the central 
star, while it is much weaker (if present) in stars with a flat disk. 
Also, a recent study of Spitzer-IRS observations of a sample of 



PDS 144 N: Our TIMMI2 spectrum of PDS 144 N clearly 
shows PAH emission at 8.6, 11.0 and 11.2 jim (see Fig.|6]). 
Comparison of the 8.6 and 11.2 fim PAH profiles with the 
different PAH classes (see Fig.^} reveals that for PDS 144 N 
the 8.6 yum PAH band belongs to class A and hence is sim- 
ilar to the ISM. In contrast, its 11.2 /um PAH profile is dis- 
tinct from both class A and B in that its peak wavelength is 
redshifted compared to both classes. Van Diedenhoven et al. 
(i2004 ) examine the 11.2 fim PAH band profile for a handful of 
HAeBe stars (other spectra were not adequate due to the pres- 
ence of crystalline silicates at this wavelength or insufficient 
SNR). All but one of these stars belong to class A for all of 
their PAH bands. More recendy, Sloan et al. ( l2005l |2007] | re- 
ported class B characteristics for four HAeBe and a T Tauri star. 
In addition, these authors found a range of peak positions for 
the 1 1 .2 fim band and relate this to the effective temperature 
of the star, i.e. the central wavelength shifts to longer wave- 
lengths as the effective temperature decreases. Tielens (I2008I I 
pointed out that this relationship can also be interpreted as a 
step function instead of a continuous dependence. In this re- 
gard, the deviation from class A seems to be connected to 
the presence of circumstellar material (Van Kerckhoven I2002t 
Boersma et al. 20081 ). However, the underlying cause of these 
variations in the 1 1.2 //m band is not well understood (see e.g. 
van Diedenhoven et al. |2004| l. With a 'central wavelength' (the 
wavelength with half the emitted flux to either side, as defined 
by Sloan et al. 120071 ) of 1 1 .289 /um and an effective temperature 
of 8750 K, this star supports this relationship (note that due the 
asymmetry of the profile, the 'central wavelength' does not cor- 
respond to the peak position). Given the class A 8.6 /im PAH 
emission, it seems that changes happened to the present PAH 
family in this source compared to that of the ISM and that these 
changes only seem to influence the 11.2 fj.m and not the 8.6 fim 
PAH bands. Observations of the CC modes (6.2 and "7.7" jum 
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Fig. 9. A comparison of the 6.2, "7.7", 8.6 and 11.2 jum PAH band profiles of HD 34700 with the PAH classes as determined by 
Peeters et al. (2002) and van Diedenhoven et al. ( |2004| l. Observations of HD 34700 are shown in thin black (both Spitzer-IRS 
and TIMMI2), the PAH classes are shown in grey shades. Clearly, HD 34700 shows unique PAH band profiles. Errors are given 
in the caption of Fig.|2l 



PAH bands, which are not observable from the ground) and a 
study of PAH emission in a larger sample of YSOs is needed to 
clarify this issue. 

No N-band silicate emission is seen for the northern compo- 
nent, which can be explained if the silicate grains located in the 
inner parts of the disk are either ( 1 ) too cold, resp. the disk lacks 
a hot optically thin surface, or (2) too large to emit at 10 //m, as 
shown by Meeus et al. (I2002I I. Our spectrum with the promi- 
nent PAH features can be explained by a flared disk model (e.g. 
Acke & van den Ancker l2004l) . which is also clearly recognised 
in the resolved images by Perrin et al. (120061) . who show that 
the disk is optically thick in the NIR. 



HD 34700: A comparison between our TIMM12 observa- 
tion and a Spitzer-IRS observation of HD 34700 is shown in 
Fig.[8l and clearly illustrates the presence of the 6.2, "7.7" and 
11.2 fim PAH emission bands. A hint of the 8.6 and 12.7 jjm 
PAH bands is also present in the Spitzer-IRS observations 
while the weaker 12.0, 13.2 and 14.2 fim bands (e.g. Hony et al. 
12001 1 and a possible 15-20 yum PAH plateau/complex (e.g. Van 
Kerckhoven et al. .2000. Peeters et al. i2004bi) are not observed 
in either spectrum. 

A detailed comparison of the PAH band profiles is shown in 
Fig.m The PAH bands exhibited by the TTS HD 34700 are 
clearly distinct. In particular, its 6.2 yum band peaks at a sim- 
ilar wavelength as that of the representative of class B, the 
Red Rectangle, but also shows excess emission shortwards of 
6. 1 yum and the long-wavelength tail is much weaker A similar 
profile is observed towards the HAeBe star HD 141569 (Sloan 
et al. 120051) . both sources showing a very strong 6.0 jum PAH 
band (6.0/6.2 band ratio of -0.33 and 0.42 for HD 141569 and 
HD 34700, respectively). Its "7.7" jum complex peaks around 
8.05 jum, redshifted compared to the Red Rectangle but, as 
noted above, class B shows a large variety and other evolved 
stars such as, e.g. IRAS 17347-3139, peak close to 8 jum. 



However, its profile is very broad, starting at similar wave- 
lengths compared to that of the Red Rectangle, but extend- 
ing up to 9 jum, thus it is more similar but blue-shifted to the 
class C profile. The latter class shows very weak, if any, 8.6 
jum emission which is also very weak in HD 34700. Similar 
broad but unique PAH profiles are also seen around the HAeBe 
stars HD 141569 and HD 135344 (Sloan et al. l2005b . suggest- 
ing a smooth transition in profiles between class B and C. The 
11.2 fj.m PAH band in HD 34700 peaks at a longer wavelength 
and also shows excess emission longwards, when compared 
to the profile of the Red Rectangle. Note that the 11.0 fim 
PAH band, usually attributed to ionised PAHs (Hudgins & 
Allamandola 119991 Hony et al. 2001), is not observed in 
HD 34700. The TIMMI2 and Spitzer-IRS profiles are rather 
similar with only small differences at the blue wing of the 
1 1.2 fim band. The 'central wavelengths' (the wavelength with 
half the emitted flux on either side, as defined by Sloan et al. 
120071) are 8.06, 1 1.35 and 1 1.37 jum for respectively the 7.7 fim 
complex, the 11.2 fim band in the IRS-SH observations and the 
1 1.2 jum band in the T1MM12 observations. While this is con- 
sistent with the results of Sloan et al. (120071 1 for the 11.2 /urn 
band, the 'central wavelength' of the 7.7 fj.m complex is con- 
siderably shorter for an effective temperature of 6000 K. 

Summarising, the PAH band profiles of the T Tauri star 
HD 34700 are clearly unique and show that PAH processing 
has occurred in this source. This is similar to PAH observations 
towards isolated HAeBe stars and a few other T Tauri stars and 
suggests that similar PAH processing occurs in protoplanetary 
disks around low-mass and intermediate-mass stars. 

Just like for PDS 144 N, the absence of silicate emission in 
HD 34700 can be explained if the silicate grains located in the 
inner parts of the disk are either too cold or too large to emit 
at 10 fim. The shallow sub-mm slope of the SED in Fig.|2]sug- 
gests grain growth for those particles located in the outer disk 
and emitting at (sub-)mm wavelengths, but this does not nec- 
essarily imply the same evolution in the more inward disk re- 



Table 4. Labels and luminosities (with errors when they were 
given in the literature) used in Figs. [T0HT2l 
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Label 


Target 


L/Lo 


Ref. 


1 


HD 72106 B 


9.9!^^,^ 


(1) 


2 


HD 98800 B 


0.58!«;il 


(2) 


3 


HD 113766 A 


4 4+0.4 
-0.4 


(3) 


4 


HD 123356 B 




(4) 


5 


HD 143006 


0.8 


(5) 


6 


HD 172555 


>16 


(6) 


7 


HD 190073 


83 


(7) 


8 


HD 319139 


n 4Q+0.06 , n qo+0.04 
"•^^-0.06 ^ "•-'-'-0.04 


(8) 


9 


CD-43 344 




(9) 


10 


FU Ori A 


466 


(10) 


11 


KK Oph A 


20!' 


(11) 


12 


MPMus 


>0.45 


(6) 


13 


PDS 144 S 


>16 


(6) 



References: (1) Folsom et al. ( |2007] |. (2) Prato et al. (|200T]l, (3) Lisse 
et al. ( 12008b . (4) neither distance nor the spectral type are certain, 
(5) Natta et al. ( I2004t . (6) lower limit when assuming a ZAMS lumi- 
nosity for the given spectral types, (7) Acke et al. ( I2005t , (8) Quast et 
al. J2000I I. (9) no luminosity class is known, (10) Green et al. (2006 ), 
(11) adapted by Carmona et al. (i2007i l. 



gions. Although much less IR excess is seen in this SED, the 
prominent PAH emission suggests that the disk is still flared 
to some extent, and not self-shadowed. Alternatively, and tak- 
ing into account that this system is a spectroscopic binary, the 
lack of warm dust could also be explained by an inner hole in 
a circumbinary disk, while the PAHs may be located further 
outwards and still emit at 10 fim. 



4. Results 

4.1. Relations between stellar, disk and silicate 
properties 

In the following we discuss the silicate shapes and study the 
silicate properties as a function of (1) stellar luminosity - be- 
tween ~0.5 and -500 L©, (2) binarity - with separations from 
approximately 0.1 to 5000 AU, and (3) infrared SED slope. To 
obtain more meaningful relations, we also consider the targets 
presented in Paper 1 and II. Apart from the three debris disks 
HD 113766, HD 123356 and HD 172555, all other targets are 
PMS objects. We will distinguish both disk types in the follow- 
ing discussion, when required. 




2.0 2.5 3.0 
Peak strength 



4.0 



Fig. 10. Silicate shape (the 11.3/9.8 ratio) in function of peak 
strength. While HD 1 1 3766 (data point ' 3 ' ) stands out due to its 
secondary generation dust, and sources of very pristine dust are 
missing in our sample, the data follow the same trend as found 
in previous studies. The target number labels are explained in 
Table |4] Typical error bars are given in the top right corner For 
comparison, data from Sicilia-Aguilar et al. (.2007i) is included 
(small diamonds). 



Silicate shape in function of peak strength: The silicate 
mass fractions (Table O show no clear dependence on stellar 
parameters. As a function of age, we see a decrease of amor- 
phous silicates and an increase of crystalline forsterite and en- 
statite, in agreement with common expectations, but no trend 
of grain growth with disk lifetime. In order to have the dis- 
cussion on more solid grounds, we use spectral indices as they 
are more robust, since they are model independent. The sili- 
cate shape, i.e. the ratio of the fluxes at 11.3 and 9.8 //m in 
the continuum normalised spectra, constitutes a good probe to 
study trends and correlations, and can be seen as a direct indi- 
cator for the grain size (Bouwman et al. 120011 van Boekel et al. 
l2003T l. 

The silicate shape in function of the peak strength is shown 
in Fig. [To] To make our results comparable to previous stud- 
ies, we follow the same definition of the silicate peak strength 
and shape as given in Sect. 4.2 of van Boekel et al. (I20051 I. 
In a weighted linear fit with relation y - a + bx, where x is 
the feature strength and y the flux ratio, values for b between 
-0.28 <b< -0.16 have been reported (van Boekel et al. l2003l 
Przygodda et al. l2003l Kessler-Silacci et al. l2007] l. We cannot 
give a mathematical relation for our sample, as it is lacking 
sources of very pristine dust (which would populate the middle 
bottom in Fig.fTOli. In order to compare with a larger sample, 
we added data from Sicilia-Aguilar et al. (I2007I I. The Vega- 
type HD 113766 (label number 3) has an isolated position in 
our plot, due to its second generation dust and its double peak 
in the spectrum. Besides this object, our sample of PMS disks 
agrees qualitatively well with the correlations found in the cited 
works: a decrease of the peak strength correlates with increas- 
ing dust processing (i.e. with a larger 1 1.3/9.8 ratio). 
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Grain processing as a function of stellar luminosity: No 

correlation of the crystalline mass fractions (Table[3j with stel- 
lar luminosity is found in our sample, nor for the derived grain 
size. Kessler-Silacci et al. (2007) have shown a correlation of 
the silicate shape with stellar luminosity, in form of an observed 
larger grain size (as witnessed by the 1 1.3/9.8 ratio) with lower 
luminosity. This trend can be understood when considering that 
the 10 yum silicate emission region lies further inward for stars 
with lower luminosity, assuming that the grain size distribution 
with radius is the same in all disks, and grain growth is most ef- 
ficient in the inner, denser regions. Our sample is much smaller 
than that presented by Kessler-Silacci et al. (2007), thus we 
cannot give a mathematical relation due to the scatter of our 
data points, but Fig. [TT] suggests a similar trend with luminos- 
ity. 



Influence of companions: We see no correlation of the de- 
gree of dust processing with stellar binarity in our sample, i.e. 
whether the time spans of dust processing and disk life time 
in binary systems would differ from those of isolated stars of 
the same age. 75% of our targets are binaries or multiple sys- 
tems. These include systems without N-band silicate emission 
(e.g. HD 3003, HD 80951, HD 155555, HD 181296), objects 
with mainly unprocessed dust (FU Ori) or large dust grains 
(KK Oph, HD 319139), and sources with all types of silicates 
present (especially HD 123356). 

For most binaries in our sample the separation (cf. Table [U is 
around 200 AU - apart from PDS 144 with -5000 AU, while 
dust in very different processing stages is found. In case of the 
three spectroscopic binary (SB) systems, silicate dust is either 
seen in small amounts or is absent, which can be explained 
by the removal of dust inside the circumbinary disk. Recently, 
Pascucci et al. ( I2008I I have shown that companions at separa- 
tions of tens of AU have no impact on the dust evolution. It 
would be interesting to study this also for a well defined target 
sample with closer companions. 



Silicate strength in function of SED slope: With ongoing 
grain growth and subsequent settling, an initially flared PMS 
disk is thought to evolve into a flat disk (Dullemond & Dominik 
I2004I I. accompanied by a change in the slope of the SED: the 
initially positive 2-25 jjm slope will decrease. During the evo- 
lution into a debris disk, the IR emission from the inner disk 
decreases when the disk is cleared from inside out. This will re- 
sult in a steeper negative SED slope with progressing evolution 
into a debris disk. When the dust has nearly completely dis- 
appeared, the slope will approach that of a blackbody with the 
stellar temperature. A dependence of silicate emission on inner 
disk evolution would be recognisable as a correlation between 
silicate strength and SED slope. There is a tendency that more 
flared (and therefore redder) disks have stronger silicate fea- 
tures. We used an average extinction curve (Savage & Mathis 
1 19791 1 with R = 3.1 and the A,, in Table [T] to calculate the de- 



Silicate shape vs. stellar luminosity 
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Fig. 11. Silicate shape in function of stellar luminosity. Objects 
with uncertain luminosity are excluded in this plot. Typical er- 
ror bars are shown in the top right corner The luminosities are 
Hsted in Tableg] 
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Fig. 12. Silicate peak strength versus infrared 2-25 /im SED 
slope (lilog / d log A). HD 113766 (data point '3') has a 
particular position due to its large amounts of secondary gen- 
eration dust. See Table |4] for the number labels. Typical error 
bars are plotted in the top right corner 



reddening of the SEDs. The slope between 2-25 fim is calcu- 
lated as 



02-25 - 



log F25 - log F2 

log .125 - log A2 ' 



(2) 



where F is the flux in Jansky at the indicated wavelengths. 
Sicilia-Aguilar et al. (I2007I I noticed a very weak to absent cor- 
relation between silicate strength and SED slope, with which 
our data agree (see Fig.fT2li. Vega-type and pre-main sequence 
stars can be clearly distinguished by their different slopes, but 
a convincing dependence on the silicate strength is not seen. 
CD-43 344 has an unusual position for a PMS source in this 
plot, but as noted earlier, its classification is not very firm. 
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4.2. Discussion 

We have presented new N-band spectra of 15 PMS and Vega- 
type stars. Besides FU Ori, these spectra were either previously 
unpubUshed or not yet available in a good quality to permit 
the analysis of the PAH bands or a dust decomposition of the 
silicate emission. HD 155555 AB is the only T Tauri star in our 
sample without an IR excess, while four of the five Vega-like 
stars only show photospheric emission at 10 //m, with excess 
emission only at longer wavelengths. For the remaining targets 
we analysed their silicate and PAH features: 

- Silicate dust: Unprocessed dust dominates the circumstel- 
lar matter around FU Ori, HD 123356, HD 143006 and 
CD-43 344. HD 123356 is considered a Vega-type ob- 
ject, and we show for the first time that the dust emis- 
sion arises from the secondary. No firm age estimation 
is available and we do not rule out that this system may 
be younger than currently believed. The N-band profile of 
HD 123356 B is rather unusual, dominated by small amor- 
phous silicates with all types of processed dust present as 
well. High-resolution NIR and MIR studies are required to 
learn more about this system. Similarly, the classification 
of CD-43 344 as a T Tauri star is not unambiguous. No dis- 
tance, luminosity class or age estimation is available. Its 
dust properties could also be explained by an evolved Mira- 
type star. Optical spectroscopy is required to determine the 
stellar class. 

Large amorphous silicates are the main dust compo- 
nent towards HD 190073, HD 319139, KK Oph, and 
PDS 144 S. No source in our sample shows highly pro- 
cessed dust, i.e. a spectrum dominated by crystalline sil- 
icates, as we found for some targets presented in Paper 1 
and II (e.g. the Herbig Ae star HD 72106 B or the Vega- 
type HD 113766 A). In agreement with previous studies, 
we find no correlation of dust processing with stellar age in 
our samples. Although HD 181296 and HD 172555 have 
a similar age and spectral type, the first source has no ex- 
cess emission in the N-band, while the last, presented in 
Paper II, shows dust in an early stage of processing. This 
could be explained by a cleared inner disk region of differ- 
ent sizes, with the disk region sensitive to N-band observa- 
tions (<10 AU) already cleared in HD 181296. 

We see no correlation of the dust evolution with binarity in 
our sample, i.e. whether the dust processing in binary sys- 
tems is retarded or accelerated when compared to isolated 
stars. 

- PAH emission: PAH bands are observed towards the 
HAeBe star PDS 144 N and the TTS HD 34700. Their pres- 
ence is often observed in HAeBe stars and commonly ex- 
plained with the excitation of PAH molecules by UV pho- 
tons. However, Li & Draine (I2002I I have shown that PAH 
molecules around cool stars can also be excited. Laboratory 
experiments by Mattioda et al. (120051) revealed electronic 
absorption bands of PAH ions up to the NIR. Thus, PAH 
emission can also be found in UV-poor environments. In 



5-35 yum spectra of a sample of 38 TTS, Geers et al. (I2006I I 
found PAHs for 3 targets (with possible further tentative de- 
tections). In their sample, T Cha of spectral type G8 repre- 
sents the lowest mass source with PAH emission. Furlan et 
al. (I2OO6I 1 found 2 cool stars with PAH emission in 5-36 fim 
spectra of 1 1 1 targets in the Taurus-Auriga star-forming re- 
gion. They show UX Tau A as the coolest TTS with PAH 
up to now, having a spectral type K5, but the PAH feature 
in this source is rather marginal. The very latest-type PAH 
emitting TTS presently is IRS 48 (Geers et al. 2007a), while 
these authors point out that its apparent spectral type (MO) 
has been highly variable in the last decade, possibly due to 
a FUOR-type accretion event. 

The PAH profiles of HD 34700 show that PAH processing 
occurs similarly to PAHs observed towards isolated HAeBe 
stars, suggesting that similar PAH processing exists in the 
protoplanetary disks around low-mass and intermediate- 
mass stars. The case is less clear for the HAeBe star 
PDS 144 N. Although the 8.6 yum PAH band profile is iden- 
tical to that in the ISM, the 11.2 /im PAH band is clearly 
distinct. 

PDS 144 constitutes a very interesting system, with prominent 
PAH bands towards the Northern source, while the Southern 
component is dominated by silicate emission. Although it is 
commonly beUeved that both objects form a binary, this re- 
mains to be proven. Perrin et al. (private communication) re- 
cently obtained optical spectra to determine more accurate 
spectral types and an improved distance estimate, to give a 
more solid answer to this question. 



5. Conclusions 

We presented new N-band spectroscopy and modelled the 
10 fj.m emission features, in order to characterise the circum- 
stellar dust in our sample. We also discovered previously un- 
known circumstellar disks. Our main conclusions are: 

1. We characterised the PAH emission in PDS 144 N and 
HD 34700. PDS 144 N exhibits mixed PAH class charac- 
teristics, while HD 34700 clearly belongs to class B. The 
latter points towards processing of PAHs, suggesting a sim- 
ilar PAH evolution in the disks of Herbig and T Tauri stars. 

2. The emission spectrum of FU Ori resembles that of 
BBW 76 in Paper I, while the FUORs V346 Nor, V883 Ori 
and Z CMa of Paper I show silicate absorption features. 
This could be an effect of the disk orientation (pole-on vs. 
face-on) or represent the two FUOR categories defined by 
Quanz et al. (.20071 ). 

3. We showed the first N-band spectrum of CD-43 344, a star 
that hosts only pristine silicates and was poorly studied. 
However, its luminosity and stellar class are rather uncer- 
tain. 

4. Silicate dust with mixed properties (sizes and crystallinity) 
is found towards HD 123356 B, while high-resolution 
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imaging is required for a further understanding of this par- 
ticular system. 

5. Within our sensitivities no N-band silicate emission is seen 
for the T Tauri star HD 155555 AB and the Vega-type ob- 
jects T] Corvi, HD 3003, HD 80951 and HD 181296. 

6. We failed to detect HD 158153 and SAO 185668 despite of 
their alleged high IRAS flux and sufficient integration time. 
We attribute this to other sources that could be present in 
the IRAS beam, causing erroneous fluxes for our target. In 
fact, there are quite a lot of other stars to which an IRAS 
flux had been attributed in the 1990s, but which later got 
revised, when higher spatial resolution was possible in the 
IR. 

7. We find no dependence of dust processing with age or 
stellar binarity, for a companion separation of typically 
>100AU. 
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Appendix A: Comments on individual sources 

A.1. FUOri stars 

FU Orionis objects (FUORs) show a variability which is char- 
acterised by dramatic outbursts in optical light, followed by a 
fading phase which can last decades. The origin of this phe- 
nomenon is not clear but it is often associated with enhanced 
accretion in low-mass pre-main sequence stars (Hartmann & 
Kenyon[l996ll. 



FU Ori is the prototype of the FU Orionis phenomenon. 
While the SIMBAD database lists a spectral type of GSIavar, 
Kenyon et al. (I2000I I changed it to GOII. In contrast to Hanner 
et al. ( I1998I I. who find an increasing MIR flux shortward of 
10 fim, the spectra of Schegerer et al. ( 120061) and Quanz et al. 
( I2OO6I 1 agree very well with ours, indicating that this object 
was not variable in the MIR between Dec. 2002 and Jan. 2006 
(or at least that a possible period with variability was not ob- 
served). No Kurucz model was plotted in the SED of FU Ori 
as explained in Sect. l3.ll but it is recognisable that the total 
IR excess of this source is larger than for all other targets in 
this work. The analysis of the silicate emission feature shows 
that it is dominated by small amorphous silicates. Beside some 
forsterite, no other dust species is found for this target in rele- 
vant amounts. 

When compared to the four FUOR spectra shown in Paper I, 
we see that FU Ori agrees well with the spectral profile of 
BBW 76, while the other 3 FUORs V346 Nor, V883 Ori and 
Z CMa show a very different feature, as it is in absorption. 
The appearance of FUOR MIR spectra could be the result of 
a geometrical effect: pole-on (emission feature) versus edge-on 
(absorption feature) view of the disk. However, given the large 
number of FUORs with silicate absorption features, this differ- 
ence cannot be explained by disk orientation alone. Therefore, 
Quanz et al. ilQQli defined two categories of FUORs: (1) ob- 
jects where the silicate feature is in absorption are likely still 
embedded in a circumstellar envelope, while (2) objects where 
the silicate band is in emission are thought to have largely shed 
their envelope, and the feature arises in the surface layer of their 
accretion disks. 



A.2. Herbig Ae/Be objects 

HD 190073 has been discussed in the scientific literature since 
more than 70 years with various classifications. Being previ- 
ously considered a post-main sequence giant, van den Ancker 
et al. ( (19981 ) describe HD 190073 as a Herbig Ae/Be star, 
what is now generally accepted. Considering its large distance, 
for which only a lower limit is known due to an imprecise 
Hipparcos parallax (van den Ancker et al. 119981 ), this unusual 
object must be much more luminous than other Herbig stars of 
similar spectral type. The optical spectrum is also quite pecu- 
liar (Pogodin et al. 12005]) . 



The SED shows an IR excess which begins in the NIR and rises 
to the far-IR. The decomposition of the N-band silicate feature 
shows a strong dominance of large amorphous silicates. Silica 
and enstatite are present in much smaller amounts, while no 
forsterite is found. 



KK Oph is a PMS binary with a 1.6" separation (e.g. Leinert 
et al. 1997 ) that shows brightness variations of the UX Orionis 
type (reddens when becoming fainter). Very controversial spec- 
tral types have been determined (between Be and FO) and the 
system was repeatedly suggested to constitute a Herbig Ae star 
with a K-type classical T Tauri companion (e.g. Leinert et al. 
[T997I Hernandez et al. [20041 Herbig l2005l l. Recently, Carmona 
et al. (12007) derived a spectral type A6Ve for the primary and 
G6Ve for the companion. Various estimates put this target at a 
distance between 160 and 310 pc, with the lower limit being 
most frequently cited. 

The N-band spectrum discussed here was already shown in 
Leinert et al. (2004), but the silicate feature was not yet anal- 
ysed in terms of dust composition. In the TIMMI2 images we 
see one isolated source, most likely just the HAe star itself, 
but we note that the MIR sensitivities were not good during 
that night. The ESO/ST-ECF Science Archive hosts another 
TIMMI2 spectrum of KK Oph, which is described in Geers et 
al. ( I2OO5I 1. These two spectra differ by about 10% in the 8 fim 
flux, while being in agreement longward of 10 fim. Both spec- 
tra, as well as the standard stars, were observed at an airmass 
close to 1.0. The distance in time between object and closest 
standard star was approximately one hour in both cases. We se- 
lected the first mentioned spectrum, since it matches closer the 
TIMMI2 Nl photometry from the Science Archive. 

The decomposition into dust components shows quite some di- 
versity. Large amorphous silicates are dominating, but silica 
and forsterite are also found in significant amounts, while en- 
statite is not seen. 

PDS 144 N: Torres et al. d 19951 1 discovered during the Pico 
dos Dias Survey the binarity of this Herbig Ae/Be star with 
a 5" separation. Spectral types of both components, given in 
Table [T] were determined by Vieira et al. (I2003I I. Recently, 
Perrin et al. (120061 ) presented comprehensive multi-wavelength 
observations and discussed the possible distance to PDS 144, 
proposing a value of 1000 + 200 pc, while previous estimates 
gave d - 140-2000 pc. Although it is commonly believed that 
both stars form a binary system, this has not yet been explicitly 
proven. The dominance of PAHs in PDS 144 N can be recog- 
nised from the NIR images in Perrin et al. (12006 ). showing an 
optically thick, edge-on flared disk, where the emission in the 
3.3 fim band clearly stands out in intensity. Until now, N-band 
spectra of both components have not yet been published. 

PDS 144 S: Resolved circumstellar emission from the more 
luminous component PDS 144 S was not detected in our and 
previous N-band images, but has recently been found by Perrin 
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et al. (private communication) at 18 /zm (0.57" diameter). In the 
N-band spectrum we see silicate features of mainly large amor- 
phous silicates, in addition to small amorphous grains and sil- 
ica. Forsterite and enstatite are absent, pointing to unprocessed 
dust, while small amounts of PAHs might be present as well. 
We note that our spectrum was taken during a night with strong 
variations in the MIR transparency, but great care was taken 
during its calibration by applying a close standard star at simi- 
lar airmass observed immediately after the science exposure. 

Both components of PDS 144 have a similar spectral type and, 
assuming they were formed together, a same age, but show 
different emission features: PAH versus silicates. The pres- 
ence/absence of PAH bands could be explained by the disk 
geometry, as we would observe much less PAH emission in 
a flat, self-shadowed disk, which prevents the UV photons 
from reaching the PAH molecules (e.g. Dullemond & Dominik 
[2004). On the other hand, we do not want to exclude the triv- 
ial explanation that the amounts of PAH in both disks could be 
different, despite the similarity of both sources. 



A. 3. T Tauri targets 

HD 34700 constitutes a quadruple T Tauri system. For the in- 
ner, short-period (P~23.5 days) spectroscopic binary, Torres 
(120041 1 had performed an orbital solution and showed that both 
stars are of similar mass and luminosity. Sterzik et al. (I2005I I 
found two more, faint, late-type components of PMS nature. 
Despite several indicators for stellar youth (IR excess. X-ray 
emission. Ha emission) the precise age has not yet been estab- 
lished, since the Hipparcos parallax is very uncertain (0.86 + 
1.84 mas). 

Sub-mm observations of HD 34700 were reported by, e.g., Syl- 
vester et al. (120011 ) and Sheret et al. (I2004I I. An N-band spec- 
trum was shown in Sylvester & Mannings (2000), but its low 
SNR does not permit a detailed study. Our T1MM12 spectrum 
is also of higher quality flian die ISO-SWS and ISO-PHOT 
data of Acke & van den Ancker (2004). As can be seen from 
the Kurucz model, excess emission starts around 5 yum and be- 
comes dominant at MIR to mm wavelengths. A PAH band fea- 
ture at 3.3 fim was shown by Smith et al. (I2004I I. 

HD 143006: As for several other PMS objects also this star 
was previously erroneously classified as pre-planetary neb- 
ula, but now it is commonly considered as a TTS (see, e.g., 
Mamajek et al. I2004I I. Contrary to the assumption of a disk- 
like distribution of the circumstellar dust. Dent et al. (i2005T l 
found single-peaked, narrow CO lines, which are more consis- 
tent with emission from an envelope than coming from a disk. 

Our SED shows a large NIR excess attributable to hot dust 
emission. Such a large NIR excess is often associated with po- 
larised emission from the circumstellar matter (Yudin [2000), 
but Hales et al. (I2006I I did not find any polarised emission and 
conclude that the disk may be too small to be resolved in their 
UKIRT observations. 



The decomposition of the N-band spectrum reveals dust at a 
very early stage of processing. Small amorphous silicates dom- 
inate the N-band feature, next to smaller amounts of silica, en- 
statite and forsterite. Moderate PAH emission is seen on top 
of the silicate feature, causing small peaks in the spectrum be- 
tween 8-9 jjm and at 1 1.2 fim. 



HD 155555 (alias V824 Arae) is a member of the /3 Pic mov- 
ing group (Torres et al. 2006) and consists of a G5IV star in a 
short-period orbit (P ~ 1 .68 days) with a KOIV-V spectroscopic 
companion (e.g. Strassmeier & Rice l2000] l. Another M4.5 vi- 
sual companion is located 33" away. The system was consid- 
ered to be a RS CVn type, however, Pasquini et al. (1991) pro- 
posed it to be a TTS. Both classifications are still being used. 
A single source is seen in our TIMMI2 images, while we were 
pointing at the position of the A and B component. The C com- 
ponent was just outside the field of view, which explains why 
IRAS obtained a higher 12 yum flux than observed in our data. 
Within our sample of T Tauri stars, HD 155555 AB constitutes 
the only target without detected excess emission above the stel- 
lar photosphere. 



HD 319139: The spectroscopic binary TTS and member of 
the /? Pic moving group (Torres et al. E506] l HD 319139 (alias 
V4046 Sgr) is known to be an emission-line star since the year 
1938. For a summary of its "early history" see, e.g., Byrne 
(119861 ). Hutchinson et al. d 19901 ) performed photometry from 
the UV to the mid-IR Q-band. From the absence of a NIR ex- 
cess, Jensen & Mathieu (1997 ) concluded there is an inner hole 
in its circumbinary disk. Quast et al. (120001 ) refined the orbital 
period to 2.421 days and determined the spectral types of both 
components to K5Ve and K7Ve. All before mentioned authors 
give distance estimates ranging between 42-83 pc. We refer in 
this work to 83 pc (Quast et al. 120001 ). as this is the most recent 
estimate. 

Our N11.9 flux (central wavelength 11.6 yum) of HD 319139 
lies about 20% below the corresponding IRAS flux, but it 
matches very well the 1 1.5 fim flux determined by Hutchinson 
et al. ( il990) . We see a rising dust continuum, on top of which 
large amorphous silicates dominate the emission feature, fol- 
lowed by small amorphous grains and some forsterite and sil- 
ica. The recently published N-band spectrum by Honda et al. 
(120061) . part of a sample of 30 low-mass stars observed with 
COMICS at the Subaru Telescope, agrees with our data. 



CD-43 344: Mannings & Barlow (fT998] ) proposed the F6V 
star HD 7151 as a debris disk candidate, by cross-correlating 
the "Michigan Catalog of Two-dimensional Spectral Types for 
the HD Stars" with the "IRAS Faint Source Survey Catalog", 
associating IRAS FO 1089-4257 with this star. Sylvester & 
Mannings (120001 ). however, showed that the more closely lo- 
cated M2 star CPD-43 142, primarily known as CD-43 344, 
can better explain the emission of this IRAS FSC source. 
Despite a significant MIR excess, this target was not yet sub- 
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ject of any further studies. Here we consider this target to be a 
TTS, based on its spectral type and pristine status of dust pro- 
cessing, but this classification is not a unique solution as, e.g., 
an evolved Mira-type star could show the same dust emission. 
Currently the distance to CD-43 344 and its luminosity class 
are undetermined. 

The 2MASS fluxes in the SED might be influenced by other 
sources, since the attributed NIR excess seems extraordinarily 
high. In the MIR, however, the target appears isolated. The de- 
composition of the silicate emission feature is dominated by 
small amorphous silicates with a less prominent contribution 
of larger grains. We did not find crystalline dust. 

A.4. Vega-type stars 

HD 3003 is a proposed member of the Tucana/Horologium 
Association (Zuckerman & Song 120041 ). In the "Catalog of 
Components of Double & Multiple stars" (Dommanget & Nys 
I2002j it is listed as a binary system with 0.1" separation and 
almost equal visual brightness of the components. Since this 
measurement is listed as being obtained in 1925, we are cau- 
tious about this information. In our T1MM12 data (with ~0.2" 
pixel resolution) only one isolated point-like emission source 
is seen. Song et al. (2001 ) give for this star an age estimation of 
50 Myr with an upper limit of 247 Myr The SED for this tar- 
get shows no IR excess emission up to the N-band, with little 
excess seen at 25 //m. 



HD 80951 constitutes a quadruple system (Dommanget & 
Nys 120021) . Oudmaijer et al. (119921 ), as well as Mannings & 
Barlow (119981) characterised it as a system with excess emis- 
sion in the IRAS passbands, but no further indication of cir- 
cumstellar matter has been described in the literature. Our 
TIMMI2 N-band spectrum is consistent with photospheric 
emission. From the SED we see an IR excess starting around 
25 fim. 

HD 109085: Stencel & Backman ([1991] ) and later Mannings 
& Barlow ( [19981 ) characterised HD 109085 (alias rj Corvi) as 
a Vega-type object. At 850 jum, a disk is found and resolved 
with a size of -100 AU (Wyatt et al. .2005) . As shown by the 
SED of this target, excess emission is insignificant at 10 fim and 
starts longward of ~17 yum, extending into the sub-mm regime. 
Our TIMMI2 N-band spectrum appears stellar. In Spitzer IRS 
data, Chen et al. (I2006I I recently detected a faint silicate feature 
around 10 fim. 



HD 123356: Not much is known about HD 123356. In "The 
Washington Visual Double Star Catalog" (Worley & Douglass 
119971 ) it is listed as having a companion 2.2 mag fainter in V- 
band, separated by 2.5" at a position angle of 126°. Mannings 
& Barlow (1998) identified HD 123356 as a candidate main se- 
quence debris disk object, by cross-correlating the "Michigan 



Catalog of Two-dimensional Spectral Types for the HD Stars" 
with the "IRAS Faint Source Survey Catalog". Sylvester et al. 
(120001) show a moderate-SNR N-band spectrum which looks 
quite different from ours. The authors, however, clarify that 
their CGS3 beam of 5.5" was pointing at the position of the 
primary (putting the secondary at the edge of the beam) and 
that not all the IR signal is detected if a substantial part of the 
emission arises from the secondary. No significant CO emis- 
sion (Dent et al. 120051 ) and only a very low degree of NIR po- 
larimetric emission (Hales et al. 120061 ) was found. 

Compared to the other emission feature sources, the TIMMI2 
N-band spectrum of HD 123356 appears atypical, since the flux 
remains rather constant between 8 and 1 1 yum. However, our N- 
band spectra were observed on three occasions (March 2003, 
February 2005 and January 2006), always showing the same 
spectral shape. From modelling the circumstellar emission, the 
dust seems to be at a rather early processing stage (or mainly 
unprocessed). As shown in Fig.|5] the emission is dominated 
by small amorphous silicates, in addition to small amounts of 
forsterite, silica, enstatite and PAH. 

While the currently available data of HD 123356 are consis- 
tent with a Vega-type star and circumstellar dust in a retarded 
stage of processing, the real nature of this system can only be 
revealed with high-resolution NIR- and MIR-observations. 

HD 181296: Mannings & Barlow ( ^9981 ) identified 
HD 181296 (ahas HR 7329, or 77 Tel) as a Vega-like star 
A brown dwarf companion at a projected distance of 200 AU 
was found by Lowrance et al. (2000). Also this target is part of 
the /3 Pic moving group (Torres et al. 120061) . 

From our TIMMI2 Nl photometry at central wavelength 
8.6 jum, the reported 12 fim IRAS flux appears overestimated, 
while our photometry and N-band spectral profile are in excel- 
lent agreement with the fluxes Mamajek et al. (120041 ) reported 
for 10.3 jum and 11.6 jum. 



A.5. Non-detections 

Two objects for which a significant MIR flux was reported 
by IRAS remain undetected in our observations, although the 
stated fluxes should have been easily observable. Most likely 
this was caused by other sources in the large IRAS beam. We 
report these non-detections to clarify the situation for these tar- 
gets in the literature. 



HD 158153: In the SIMBAD database an IRAS flux of 3.83 Jy 
at 12 jum is given. This source was also mentioned in the 
work by Oudmaijer et al. ( 119921 ). With a lOcr sensitivity of 
25 mJy/hour we did not detect a source at this position within 3 
minutes of integration time, suggesting a flux limit <1 10 mJy. 
Several IRAS sources are located nearby, of which the 11" 
close object IRAS 17251-2320 seems to emit the largest flux. 
We find a rising dust continuum for that source, but do not 
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present it here in detail, as its nature (stellar or not) remains 
unclear. 



SAO 185668 is listed in SIMBAD with an IRAS flux of 
1.61 Jy at 12 fim. Malfait et al. ( 119981 ) considered it a Vega- 
like star based on its IRAS photometry. With a poor lOcr 
MIR sensitivity around 100-150 mJy/hour no signal was found 
within 12 minutes on-source time, resulting in a flux limit about 
<0.3 Jy. 



